Scattering and bound states of two polaritons in an array of coupled cavities 
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We develop an analytical approach for calculating the scattering and bound states of two po- 
laritons in a one-dimensional (ID) infinite array of coupled cavities, with each cavity coupled to 
a two- level system (TLS). In particular, we find that in such a system a contact interaction be- 
tween two polaritons is induced by the nonlinearity of the Jaynes-Cummigs Hamiltonian. Using 
our approach we solve the two-polariton problem with zero center-of-mass momentum, and find 
ID resonances. Our results are relevant to the transport of two polaritons, and are helpful for the 
investigation of many-body physics in a dilute gas of polaritons in a ID cavity array. 

PACS numbers: 



I. INTRODUCTION 

In the recent years, the investigation of the physics 
in one-dimensional (ID) [H-flsj and two-dimensional 
(2D) p^ - [23j array of coupled cavities has attracted a 
lot of attention. It is predicted that, such systems can 
be used in quantum information processing [1,-5] as well 
as the quantum simulation of many-body systems, e.g., 
the quantum phase simulation [l0|-|20|, quantum Hall 



effect [2l| and Bose-Einstein condensate For the 



few-body physics of coupled cavities, many authors have 
studied the single-photon transmission in a ID cavity ar- 
ray coupled with a single atom and the dynamics of 
a single polariton in a ID cavity array with each cavity 
coupled to an atom 0, 0] ■ Recently, bound states of two 
polaritons in such a system with finite number of cavities 
was studied by Wong and Law by direct numerical diag- 
onalization of the Hamiltonian [8i]. However, the scat- 
tering and bound states of two polaritons in an infinite 
array of coupled cavities remain to be investigated. 

In this paper, we develop an analytical approach for 
calculating the scattering and bound states of two po- 
laritons in an infinite array of coupled cavities with each 
cavity coupled to a two- level system (TLS), which can 
be either a natural or an artificial atom. In particular, 
we find that in such a system there is an effective contact 
interaction between two polaritons. Similarly to the pho- 
ton blockade phenomenon [23], this effective interaction 
is also due to the nonlinearity of the Jaynes-Cummigs 
(JC) Hamiltonian. With our approach, one can treat the 
problem with standard techniques of quantum scattering. 
We further derive the scattering coefficient and bound- 
state energy of two polaritons with zero center-of-mass 
momentum. The resonance phenomena induced by the 
weakly bound states are also investigated. Our results 
can be directly used in the research of the transport of 
two polaritons. It is also helpful for the investigation of 
the many-body physics in a dilute gas of polaritons in a 



ID cavity array. 

This paper is organized as follows. In Sec. II we derive 
the effective interaction between two polaritons. Based 
on this result, the analytical method for calculating two- 
polariton scattering states and bound states is developed 
in Sec. HI. In Sec. IV we solve the two-polariton problem 
with zero center-of-mass momentum with our approach 
and analyze the resonance phenomena. We conclude and 
discuss these results in Sec. V. 



II. TWO-BODY PROBLEM OF POLARITONS 
IN CAVITY ARRAY 

We consider a one-dimensional (ID) infinite array of 
coupled single-mode cavities. In each cavity there is a 
two- level system (TLS), which interacts with the photon 
in the cavity. 
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FIG. 1: (Color online) Schematic illustration of a one- 
dimensional infinite array of coupled single-mode cavities, 
with each cavity coupled to a TLS. The single-photon hopping 
intensity between the adjacent cavities is ^, and the coupling 
intensity between the TLS and the photon in the same cavity 
is g. 

In the interaction picture, the Hamiltonian of the sys- 
tem is given by 

oo CXD OO 

H = A J2 '^i^e^+a E E cii+.an+h.c. , 
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where A = — wl is the detuning between the TLS 
and the cavity mode, with uja and oj^ the frequency of 
the TLS and the photon, respectively. In Eq. (HJ ^ is the 
hopping intensity or the inter-cavity couphng strength, 
a„ and ajj are the annihilation and creation operators of 
the photons in the n-th cavity respectively, and g is the 
coupling intensity between the TLS and the photon in the 
n-th cavity. Without loss of generality, here we assume 
that g and ^ are real numbers. The Pauli operator cr,|j'' 
{iJ = e,5) is defined as a^f = with \g)["'^ and 

16)4"' the ground and excited state of the TLS in the rt-th 
cavity, respectively. For the convenience of our following 
discussion, we further define |0)i"'' as the vacuum state 
of the n-th cavity. 

We can define the "vacuum" state of our total system 

as 



n 



(2) 



with K €: [—2tt, 2tt) the total momentum of the two po- 
laritons, and the coefficients ti,di,pi describe the relative 
motion of the two polaritons, and satisfy 



pi ^p-u ti = i-;, ^0 = 0. 



(6) 



We can further define the coefficients di± = (d/ ± /2 
which satisfy di± = ±d^i±. 

It is apparent that the quantum state j^") is described 
by a 4-dimensional vector 

pi = ipi,di+,di-,tif (Z = 0,±l,±2,...), (7) 

which can be considered as the "spinor wave function" 
of the relative motion of two polaritons in the state \'^). 
Thus, the condition ^ for the coefficients can be re- 
expressed as 



Here the matrix T is defined as 



(8) 



where all the atoms are in their ground level and there 
is no photon. Both the creation of photon and the exci- 
tation of TLS can be considered as the polaritons of the 
system. We can define the number operator of polaritons 
as 



N = 



00 

E 

n— — oc 



(3) 



It is clear that [H, iV] = and the polariton number is 
conserved. In case of iV = 1, there is a single polariton in 
our system, and Hamiltonian H has eigen-states jC) 
or bI,\G) with eigen-energies eAk,Bk — (2(^cosA; + A ± 
^(2^cos/c- A)2-hV)/2. Here 

00 

4(4)= E 

n— — 00 

(4) 

are the creation operators of the excitonic polariton of 
kind A ov B with momentum k. The coefficients {k) 
and 77^ g (k) can be obtained straightforwardly from the 
eigen-equaiton of H. 

In this paper we consider the two-body problem of po- 
laritons in our system, i.e., the quantum dynamics of our 
system in the subspace with N = 2. In such a subspace, 
due to the translation symmetry of our system, the total 
momentum of the two polaritons is conserved, and thus 
the eigen-state of the Hamiltonian H can be written as 



l*> = 



m,n— — 00 
00 



m,n— — 00 
00 ^ 



m,n—— 00 



V2 



1000 
0100 
00-10 
0001 



(9) 



A straightforward calculation shows that, the 
eigenequation — E\'^) of the Hamiltonian H can 

be re-written as 



(10) 



Here for each given value of (/,/'), hI^j) and V;_/' are 
4-dimensional matrixes and defined as 



H 



(0) _ 

Ll' ~ 



1,1' 
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and 
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VO -V2g 



V2.g 


/ 



i',0, 



(11) 
(12) 



with the matrixes A, B and C given by 



A 



/ V2.g \ 

V2.g A %/2g 

A 

V2g 2A y 

2 

10 

10 





B = 2^cos(X/2) 



C = -2£,sm{K/2) 



/O 

10 

10 

Vo 



(13) 



(14) 



(15) 
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The above discussion shows that the Schrodinger equa- 
tion for the relative motion of two polaritons in our sys- 
tem is equivalent to Eq. (fTOl) and the boundary condition 
for the coefficient /3;. It is obvious that Eq. (|10l) has 
a same form as the stationary Schrodinger equation for a 
scattering problem. Then we can understand the matrix 
'H.f^i, as the "free-Hamiltonian" of the two-polariton rela- 
tive motion, and V;^;/ as the "interaction" between these 
two polaritons. Therefore, we can use the standard tech- 
nique of quantum scattering problem to solve Eq. (|10p. 
and then find the scattering states and bound states of 
the two polaritons in our system. 

Now we comment on the physical picture given by the 
inter-polariton interaction ~Vi^i' . It is clear that, the sym- 
bol I {I') is an abbreviation of m — n in Eq. ([S]), and thus 
describes the relative position of two polaritons. Since 
Vi^i' takes non-zero value only when / = /' = 0, it is a 
two-polariton contact potential. Namely, two polaritons 
interact with each other when they are in the same cav- 
ity. To understand the physical meaning of this contact 
potential, we consider a extreme case where the photon 
cannot tunnel between different cavities. In that case, 
the total Hamiltonian becomes H = J^'jL^oo -^/c ''^ith 



of and find the "free-motion" state of the two po- 



laritons. Due to the translation symmetry of H 



basic solution of Eq. ([TS]) takes the form e'^'F (q) with 
q G [— TT, tt) the relative momentum of the two polaritons 
and the ^-independent vector F (q) satisfies 



[A + Bcosq + Csmq]F{q) = EF {q) . 



(17) 



It is obvious that, for a given value of q, Eq. (IT7| has four 
solutions for F{q). We denote these solutions as Fa (q) 
with a taking the values AA, AB, BA, BB. For a physi- 
cal meaning of these solutions, see the final paragraph in 
this section. The analytical expression of Fa{q) is given 
in Appendix A. Straightforward calculation also shows 
that, the eigen-energy E{a,q) with respect to e^^'F^ (q) 
can be expressed as 



E (a, q) = £u (Ao, 6i) + £^ (Aq, 62) 



(18) 



Here the symbols u, v can take the values A, i?, and 
related with a via the relationship a = uv. In Eq. 
(IT51) we also have Aq = 2^cosq'COs(i4r/2) and 61^2 — 
A q= 2^singsin(i4r/2).The function 8A,B{x,y) is defined 
as 



the JC Hamiltonian H 



JC 



h.c. It is 



(n) 

clear that Hj^ can be diagonalized in the subspace with 



1 1 / ^ 

SA,B{x,y)^ -ix + y)±-\/{x~y) + ig^ 



+ a]aj = n {n = 0, 1,2...), and the relevant eigen- We further define the vector /S^"^ {a, q) as 



energies are s^f^{n) = (A ± A'^ + Ag'^n)/2. If two po- 
laritons appear in two different cavities, the total energy 
of the two polaritons can take the (i, j)-independent value 

^JC (1) + ^jc(l)- However, if the two polaritons appear 
in the same cavity, the energy of the two polaritons take 

the value £^^^(2) ^ £jc (1) + Ejc (!)• Therefore, due to 
the nonlinearity of the spectrum of JC Hamiltonian, the 
energy of two polaritons changes when they are in the 
same cavity. That is the origin of the effective contact 
interaction V;^;/ of two polaritons in our system. 



III. SCATTERING AND BOUND STATES OF 
TWO POLARITONS 



In the above section, we find that the two-polariton 
problem in the ID cavity array is described by Eq. (jlOp 
with boundary condition In this section we show our 
approach for solving Eq. ([TU]) and derive the scattering 
states and bound states of two polaritons. 



A. Scattering states 

Now we calculate the two-polariton scattering state. 
To this end, we first analytically solve the eigenequation 



(a,<z) = -[e'^'F„(<z) 

,(0) 



-'''TFa (g)] 



(19) 



(20) 



It is easy to prove that (a, q) satisfies both Eq. (fT6)) 
and the boundary condition ([8|), and then can be consid- 
ered as the "free-motion" state of the two polaritons. 

Now we consider the scattering wave function 
Z?!"'^'' (a, q) with respect to the incident wave function 
'(")(a,g). «W 



equation 



/3; is given by the Lippmman-Schwinger 



(a, g) {a,q) 



EgC!v,,,,/3(+) (a, 9), 



where the Green's function G^^) i, 
equation 



(21) 

is the solution of the 



I" 



[E{a,q) 



H 



(0) 

1,1" 



(22) 



with I the 4-dimensional identical matrix. The straight- 
forward calculations (see, e.g., chapter 9 of Ref. E25!|) 
with the Lippmman-Schwinger equation (j2ip show that 

l3l~^^ {a,q) takes the form 



/3|+)(a,g) = /3}°^a,g)+^/(7^a,g)e^^-'F,(A^) 



,(0) 



(0) 



(16) 



/3|+'(«,g) = r/3iV(a,q) 
("'9) = (5p,s+,0,0) 



il>0), (23) 
(/<0), (24) 
(25) 
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Here F~^{\j) is given by Eq. (jAll) and the parameters state are determined by the eigenequation 
(7, A^) satisfy E{j,Xj) — E{a,q). Namely, A-y can be 



obtained via the equation 

det\A + BcosXy+CsmX^-IE{a,q)\^0. (26) 



In addition, A^ also satisfies the conditions 



ImA^, > 0; 



(27) 



T7T-E{j, A^) > when ImA^, = 0. (28) 

(JA.y 

It is easy to prove that Eqs. (I26II28P have three solutions. 
Thus, the summation in Eq. (1231) includes three terms. 
In Eq. ([25)1 . the factor f{a ^ ct,q) is the elastic scat- 
tering coefficient, while /(7 4- a,q) for the term with 
7 7^ a and ImA^ > is the inelastic scattering coefficient 
which describes the inter-channel transition induced by 
the scattering process. 

For a given incident wave function /3|*''' {a,q), we can 

obtain the scattering wave function P^^^ (a, q) with the 
following two steps. First, solve Eqs. (1261128^ and find the 
three solutions for (7, A-y). Second, substitute expressions 

and the 



p3ll25)) into equations Ei'(Hl_ 



(0) 



with / = 0, 1, and obtain the values of s+, Sj 
values of the three coefficients /(7 a, q). 

In the end of this subsection, we discuss the physical 
meaning of the two-polariton scattering state. To this 
end, we first consider a state |$) = At,,^, ^A\^,„_ \G) 



K/2+q^'^K/2-q' 

dH). The physi- 



with the operator defined in Eq 
cal meaning of the state |$) is that there are two ex- 
citonic polaritons of type A with momentums K/2 -\- q 
and K/2 — q. It is apparent that the state |$) can 
be written in the form of Eq. ([5]). We denote the (3- 
coefficient or the relative wave function of |<I>) as /3i{^). A 
straightforward calculation shows that we have — 

Z?!"'' {AA,q). Furthermore, it can also be proved that 
£A,K/2+q +£A,K/2-q = E{AA,q). Namely, the total en- 
ergy of the two excitonic polaritons in the state \^) is 
the same as the energy of the incident state with wave 
function f^^^^ {AA,q). Thus, the incident wave function 

/3;°'' {AA, q) can be considered as the relative wave func- 
tion of the "free motion" of two excitonic polaritons of 
type A with total momentum K and relative momentum 
q. Similar discussions can also be done for the states 
with two excitonic polaritons of the other types. There- 
fore, Pi^^ {uv,q) and (3^^'^ iuv,q) with u,v — A, B can 
be considered as the relative wave function of the free 
motion and scattering state of two excitonic polaritons 
of types u and v with total momentum K and relative 
momentum q. 



Bound states 



Now we consider the bound states of two polaritons. 
The energy Ei, and the wave function /Sf''-' of the bound 



and the boundary condition 



hm /3,^^' = 



(29) 



(30) 



Similar as in the above subsection, the solution of Eqs. 
\30\i takes the form 



/Sf' = ^6,e'^^'f^„^(A,), 



(b) 



(l > 0) (31) 

{I < 0) (32) 
(33) 



where is defined in our above subsection and {aj, Xj) 
is given by the relationship E{aj,Xj) = Ei, and satisfies 
the equations 



det |yl + S cos Aj + C sin 



A, 



-IEb\ 
ImAj 



= 0; 
> 0. 



(34) 
(35) 



In addition, substituting Eqs. pm33p into equations 
Eq. (|29p with ^ = 0, 1, one can get the homogeneous lin- 
ear equations for the coefficients 61,2,3, bp and When 
there is a two-polariton bound state, the determinant of 
the coefficient matrix of these equations should be zero. 
From this condition, we can obtain the energy Ei, of the 
bound state. Substituting the value of Eh into these 
linear equations and using the normalization condition 
J2i pf'^^ pf'^ — 1, we can obtain the coefficients 61,2,3, bp 
and 6+. 



IV. TWO-POLARITON PROBLEM WITH ZERO 
CENTER-OF-MASS MOMENTUM 

In the above section, we show our analytical approach 
for the calculations of two-polariton scattering states and 
bound states in a ID cavity array. With our method, one 
can derive the scattering coefficient and bound states of 
two polaritons in the systems with any detuning A and 
coupling parameters (5,0- Now we show the results of 
our calculations. For simplicity, in this paper we only 
consider case with the total momentum K = Q. 

As discussed above, the energy of the two-polariton 
scattering state P^'^'^ {uv,q) with u,v = A, B is, E{uv,q) 
defined Eq. (fTS]) . In the following, we call the region 
of value of E{uv,q) with q G [— 7r,7r) as "ww-band". 
Then the energies of the scattering states are located in 
AA-, AB-, BA-, and i3i3-bands. Furthermore, we have 
E{AB,q) = E{BA,-q). Thus the AS-band and BA- 
band totally overlap with each other, and we have three 
energy bands for the two-polariton scattering states. In 
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FIG. 2: (Color online) The energy bands of scattering states 
(regions with shadow) and eigen-energies of bound states (red 
circles) as functions of the detuning A. The figures are plot- 
ted for the cases with total momentum A' = and photonic 
hopping intensity ^ = —0.2(7 (a), — 0.5(; (b) and = — 0.75(? (c), 
with g the photon-TLS coupling intensity. 
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FIG. 3: (Color online) Norm-square of two-polariton scatter- 
ing coefficient as a function of relative momentum q of the 
incident state with total momentum K = Q and photonic 
hopping intensity ^ = — 0.2(;. (a) Norm-square of the elas- 
tic scattering coefficient /(AA <— AA^q). (b) Norm-square 
of the elastic scattering coefficient f{BB BB,q). (c) 
Norm-square of the elastic and inelastic scattering coefficient 
f{AB -It- AB,q) and f{BA -ir- AB,q). In our cases we have 
f{AB -t- AB,q) = f{BA <~ AB,q), and then denote the 
value of the two amplitudes as Fab in the figures. 



Fig. 2, these three bands are shown as the regions with 
shadow. 

We further find that, in our system there are two bound 
states with wave functions (3^''^^ and /^^''^^ and energies 
Ebi and i?f,2 which are located in the gap between AA- 
and AB-band, and the one between AB- and SB-band, 



respectively. In Fig. 2 we demonstrate the energies E{,i 
and E{,2 as functions of the detuning A. 

In the large-detuning cases with |A| >> the 
effects of photon-TLS coupling become weak. As a re- 
sult, the creation and annihilation of photons, as well as 
the quantum transitions between the ground and excited 
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states of TLS, can be adiabatically eliminated. In addi- 
tion, the effective coupling between photons in different 
cavities, as well as the one between the excited TLSs in 
different cavities, can appear in our system. Then one of 
the two bound states can be approximated as the two- 
photon bound state, and another one becomes the bound 
state of two excitations of TLS. Furthermore, the ener- 
gies of the two bound states become very close to the 
borders of the energy bands of the scattering states. 

The above analysis are verified by our quantitative 
calculation. Our results show that, when A < and 
|A| >> |g|,|^|, the energy Eti approaches to the lower 
limit of the AA-band (Fig. 2), and we have Z?}^^^ ~ 
(p;, 0, 0, 0)"^. Therefore, the two-polariton bound state 

l^'bi) satisfies l^-bi) w Y.,n,iPi"'lnaln+i\G) , and it is ap- 
proximately a two-photon bound state. On the other 
hand, when A >> j^l, |^|, the energy Ebi approaches to 
the upper limit of the AB-hand (Fig. 2), and we have 

PI''^^ « (0, 0, 0, a/)-^. Therefore, the two-polariton bound 



state I'^bi) satisfies I^E" 



bi/ « E™,z*i<^e"V^"+'^|G), and 
becomes approximately a bound state of two excitations 
of TLS. 

Similar analysis can be done for the bound state \'^b2) 
with wave function /3|''^'' and energy £'f,2- As shown in 
Fig. 2, the energy Eb2 approaches to the upper limit of 
BB-ha.nd and the lower limit of AB-hand in the limits 
A — > ±c», respectively. Furthermore, \'^b2) is approx- 
imately a two-photon bound state when A >> Ig], |^|, 
and approximately a bound state of two excitations of 
TLS when A < and |A| » \^\. 

In the ID low-energy scattering problem of two non- 
relativistic particles in the continuous space (e.g., the 
scattering problem of a non-relativistic particle in a delta 
potential gS{x) with 5 < 0), it is well-known that, when 
the energy of the bound state is close to the lower limit 
of the scattering energies, the scattering coeflScient (re- 
flection coefficient) becomes very small. In our system, 
we find similar resonance results. As shown above, in 
our cases the band of the scattering energies has both 
the lower border and the upper border. We find that 
when the bound-state energy is close to either of these 
two borders, the relevant scattering coefficient / becomes 
very small. 

For the elastic scattering coefficient f{AA ^ AA,q) 
of the scattering state in the AA band, as shown in 
Fig. 3(a), when A < and |A| >> \g\,\^\ and £^61 is 
close to the lower bound of AA-band, we have f{AA 
AA,q) w 0. On the other hand, when A >> \g\, |^| and 
Ebi is relatively far from the lower bound of ^^-band we 
get the result |/(AA ^ AA, g)| « 1. Our results show 
that, the effective interaction between two excitonic po- 
laritons of kind A is negligible in the limit A — — 00, and 
becomes strongly repulsive when A ^ cx). 

For the elastic scattering coefficients f{BB ^ BB, q), 
f{AB <— AB,q) and the inelastic scattering coeffi- 
cient f{BA •(— AB,q), as shown in Fig. 3(b), we have 
f{BB ^ BB,q) sa when A >> |g|,|^| and Eb2 is 



close to the upper bound of _Bi?-band. Finally, for the 
scattering states in the AB{BA) band, according to our 
calculation (Fig. 2), the energy Ebi is close to the up- 
per bound of the AB-hand when A >> while 
Eb2 is close to the lower bound of the AB-hand when 
A < and |A| >> The corresponding reso- 

nance phenomenon is illustrated in Fig. 3(c), where it is 
shown that in both of these two regions we always have 
f{AB ^ AB, q) f{AB ^ BA, q) « 0. 



V. CONCLUSION AND DISCUSSION 

In this paper we derive the effective interaction be- 
tween two polaritons in a ID cavity- array coupled to 
TLSs, and provide an analytical method for the calcula- 
tion of scattering and bound states in such a system. We 
find that in a cavity array there is an effective contact in- 
teraction between two polaritons, which is induced by the 
nonlinearity of the JC Hamiltonian. The interaction is 
totally determined by the photon-TLS coupling strength 
g, and independent of the photonic hopping intensity ^. 
For two polaritons with zero center-of-mass momentum, 
we find that there are two bound states in the gaps be- 
tween the energy bands of the polaritons, and the ID res- 
onance phenomenon can appear when the photon-TLS 
detuning is large enough. Our result is helpful for the 
research of the few-body and many-body physics of po- 
laritons in cavity arraies. In particular, for the dilute po- 
lariton gas where the average polariton number in each 
cavity is much smaller than one, the many-body physics 
is dominated by the two-polariton contact interaction. 
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Appendix A: The expression of Fa (q) 

In this appendix we provide the analytical expression 
of the vector (q) (a = AA, AB, BA, BB) defined 
in Sec. III. With straightforward calculation, we solve 
Eq. p7|) analytically and get the expression of Fuv (q) 
{u,v = A, B): 



Fuv (q) 



/ r^(l)^(2) 

^gu ^gv 

1 /^/^(l)/^(2) , ^(1)^(2)^ 

1 ?/^(l)/^(2) _ ^(1)^(2)^ 
\^eu ^gv ^gu ^ev 



\ 



(Al) 
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where the coefficients c''e(g) a{b) 

(1.2) 



are defined as 



e(g)[A(B) = fe{g)A(B) {^1,2 - Aq) with 5i^2 and Ao de- 



c 

fined in Sec. IIIA. Here the functions fe(g),A(B) {x) are 
defined as 



fg,A(B) {x) 
fe,A{B) {x) 



4r 



(A2) 



(A3) 
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